Change in Order of Phase Transitions on Fractal Lattices 
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We re-examine a population model which exhibits a continuous absorbing phase transition which 
belongs to directed percolation in 1+1 dimensions and a first order transition in 2+1 dimensions and 
above. Studying the model on fractal lattices, we examine at what fractal dimension 1 < df < 2, the 
change in order occurs. As well as commenting on the order of the transitions, we produce estimates 
for the critical points and, for continuous transitions, some critical exponents. 
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Continuous nonequilibrium phase transitions continue 
to be an area of great interest (see [l|, Q for recent 
reviews). Perhaps the main reason for this is due to 
the observed universahty that many such models dis- 
play. As in equilibrium phase transitions, models be- 
longing to the same universality class have identical crit- 
ical exponents and the scaling functions become identi- 
cal close to the critical point. By far the largest uni- 
versality class of nonequilibrium phase transitions is di- 
rected percolation (DP). Indeed, it has been conjectured 
by Janssen and Grassberger that all models with a scalar 
order-parameter that exhibit a continuous phase transi- 
tion from an active state to a single absorbing state be- 
long to the class d, 0] ■ 

Most of the research on phase transitions has been 
conducted on lattices of integer dimension. Motivated 
by the apparent dependence of critical exponents on 
dimension, as well as topological features such as ram- 
ification and connectivity, work has also been carried 
out on lattices of fractal dimension d{ [^, 0|. In 
this paper, we examine a model where the order of the 
phase transition changes from continuous to first-order 
for 1 < df < 2. We investigate the model on fractal 
lattices to examine at what dimension the change in 
order occurs d[^. 

The model: In a recent paper Q, we introduced a 
general population model with the reactions 



2A + 



3A, A 



and A(f> < — > (f>A, (1) 



for an individual A. The model is simulated on d- 
dimensional square lattices of linear length L where each 
site is either occupied by a single particle (1) or is empty 
(0). A site is chosen at random. With probability pd the 
particle on an occupied site dies, leaving the site empty. 
If the particle does not die, a nearest neighbour site is 
randomly chosen. If the neighbouring site is empty the 
particle moves there, otherwise, the particle reproduces 
with probability pb producing a new particle on another 
randomly selected neighbouring site, conditional on that 



site being empty. A time-step is defined as the number 
of lattice sites N — L"^ and periodic boundary conditions 
are used. Assuming spatial homogeneity, the mean field 
(MF) equation is easily found to be 



dt 



Pb(l -Pd)p(t)'(l -pW) -PdP(t). (2) 



This has three stationary states. 



Po = 0, p± 



1 ± 



4pd 



Phil -Pd) 



(3) 



For Pd > pb/(4 + Pb), p± are imaginary resulting in 
Po being the only real stationary state. We therefore 
define the critical point pd — Pd^ = ?5b/(4 +Pb), marking 
the first-order phase transition between survival and 
extinction of the population. As we showed in our 
paper Q, from monte carlo (MC) simulations, the 
model displays a continuous phase transition in 1+1 
dimensions and belongs to DP. In higher integer di- 
mensions, consistent with the MF prediction, the model 
displays a first-order phase transition. The difference in 
order of the phase transition between the MF and MC 
simulation results in 1+1 dimensions only is likely due 
to the larger correlations present in lower dimensions. 
Here, an individual is very likely to be able to find a 
mate on a neighbouring site. In 1+1 dimensions then, 
the reproduction rate is more accurately proportional to 
p{l ~ p) rather than p^{l — p). 

For the continuous phase transitions, we examined 
the population size n{t) and survival probability P{t) 
after beginning the simulations from a single seed - two 
adjacent particles. We expect the asymptotic power law 
behaviour 191 



i{t) 6c and P{t) 6c t~ 



(4) 
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at the critical point. Such time-dependent simulations 
are known to be more efficient than steady state simu- 
lations which are also more prone to finite-size effects 
d El- Fig. [1] a) shows the resuffs for the (1+1)- 
dimensional case where the gradient in the log-log plots 
gives the DP exponents f] and S. For the first-order 
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FIG. 1: a) Power law behaviour at the continuous phase tran- 
sition in 1+1 dimensions. The hashed hnes show the DP val- 
ues rj = 0.313686 and 5 = 0.159464 [TJ. b) The values of 
Pdc (L) with the hashed line of best fit which extrapolates the 
data as L ^ oo. The inset shows the double-peaked structure 
at the critical point for a finite L. 



phase transitions, we began our simulations from a fully- 
occupied lattice and observed a double-peaked structure 
in the histogram of population density M{p). Finding 
the value of pd which equated the size of the two peaks 
at po and p+ gave the value of (-^) due to the phase 
coexistence that occurs at first-order transitions. Extrap- 
olating the data as L — )■ cx) gave a value for the critical 
point as shown in Fig. [l]b). 

Having described how we examined our model in 
integer dimensions, we turn now to describing our 
methodology for the fractal lattices. 

Methodology and results: Sierpinski carpets have 
been used widely to provide a generic model of fractals 
to study physical phenomena in fractal dimensions (see 
for example [1, El IH, H, El El)- They are formed 
by dividing a square into P sub-squares and removing 
{P — Noc) of these sub-squares from the centre This 
procedure is then iterated on the remaining subsquares 
and repeated k times. As k — > cxd a fractal structure. 



FIG. 2: Snapshots of the population on the fractal lattices 
S'C(5^16,3) (left) and S'C(3^8,4). The left fractal with 
di = log(16)/ log(5) ~ 17227 was initiated with a single seed, 
whereas the right lattice with dt = log(8)/ log(3) ~ 1.8928 
began from a fully-occupied lattice. 



denoted by SC{P , Noc), is formed. For finite k, how- 
ever, we denote the structure SC{P , Nqc, k) which has 
N = sites. The Hausdorff fractal dimension of the 
structure as k — > oo is then df — In(A'oc)/ lii(^). Using 
different values of I and Noc enables us to use lattices of 
different fractal dimension d{ where 1 < df < 2. 

Using these fractal-dimensional lattice structures, we 
carried out our simulations as before to determine the 
values for p^^ and, where a continuous phase transition 
occurred, the critical exponents. For continuous phase 
transitions, we used time-dependent simulations, begin- 
ning from a single-seed. Here, however, the position of 
our initial seed may well affect the results. Two adjacent 
particles next to a large hole would clearly be less likely 
to survive than two particles surrounded by empty sites. 
We therefore randomly picked two adjacent sites for each 
run and, as before, made an average over all runs. It is 
also clear that, unlike before, we are unable to begin the 
simulations with a single-seed at the centre of the lat- 
tice due to the way in which the lattices are constructed. 
Due to the random position of this single-seed, we are 
not able to make the lattice sufficiently large so that the 
particles never reach the boundary. As a compromise, we 
use very large lattices and periodic boundary conditions. 
For the first-order phase transitions, we examined the 
histograms of population density having started from a 
fully-occupied lattice as described earlier. Example snap- 
shots of both of these approaches are shown in Fig. [21 

There are many problems associated with examining 
critical behaviour on fractal lattices. Issues of boundary 
conditions and the initial location of the single seed have 
already been mentioned. Further, there is the obvious 
difficulty in the fact that finite lattices are not truly frac- 
tal. Indeed, for the histogram approach, we are required 
to find pd^{N) for different lattice sizes which we achieve 
by increasing k. Such an approach, however, changes 
the structure of the lattice and can therefore affect the 
results. The effects of these difficulties and inaccuracies 
are minimised here since we are primarily interested in 
the order of the phase transition at the different dimen- 
sions. 
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FIG. 3: a) Plots of n{t) and P{t) at the critical point for 
fractal dimensions 1.5573 and 1.7227. The hashed lines give 
the estimated values for the exponents outlined in Tab. U 
The inset shows R^(t) for df ~ 1.5573 only. The hashed 
line has gradient 2/z where z is obtained using the scaling 
relation (|6]) and the obtained values for t; and S. b) The 
double peaked histogram of population density indicating a 
first-order phase transition for df ~ 1.8928. The inset shows 
the predicted values of pdc (N) and an extrapolation of these 
results for N oo. 



For df ~ 1.5573 and d{ ~ 1.7227, we find continu- 
ous phase transitions shown by observing power law be- 
haviour at the critical point. Plots of both n{t) and 
P{t) are shown in Fig. [3] a). For df ~ 1.5573, we 
used S'C(9^32,4) which has just over 10*^ sites. The 
data was obtained from over 1.6 x 10^ independent runs. 
Comparing the plots for this dimension with those for 
d = 1, we observe significantly larger corrections to scal- 
ing. Whereas power law behaviour was obtained after 
~ 10^ time steps for the d — 1 case, for df ~ 1.5573, 
we have to wait ^ 10^ time steps and ~ lO'^ time steps 
for df ~ 1.7227. For this latter case, we therefore had to 
increase the number of time steps as well as the number 
of independent runs to 2.25 x 10^. The simulations took 
over three months of computer time on Imperial College 
London's HFC for each parameter value. As the 




FIG. 4: Possible power law behaviour for df ~ 1.7927. The 
inset shows the lack of the double-peaked structure. 



crossover in the order of transition is approached, these 
corrections to scaling are likely to increase further in size. 
This, then, represents a further challenge in obtaining 
accurate values for the critical points and exponents as 
df df^. 

From the other end of the interval 1 < df < 2, for df ~ 
1.8928, no power law behaviour was observed, rather, the 
double-peaked structure of the histogram of population 
density. To increase the size of our lattice, we increased 
the value of k and found the value of pd^{N) for each 
case. Extrapolating the results for infinite system size, 
we obtained an approximation for the critical point as 
shown in Fig. [3]b). 

For df ~ 1.7925, the histogram plots do not show the 
double-peaked structure close to the critical point. Plots 
of n{t) and P{t) appear to show power law behaviour, but 
with very large corrections to scaling. The plots shown in 
Fig. IHrequired over 5 x 10^ runs, taking over six months 
of computer time. On closer inspection, the values of 
77 = 0.19 and S — 0.46 that we obtained at this dimension 
appear to be unlikely. Both vales are the wrong side of 
the DP values for d — 2, indicating that the plots actually 
show supercritical behaviour. Due to the large amounts 
of time required to produce such plots, we predict that 
the phase transition is continuous for df 2± 1.7925 but 
are, now at least, unable to provide strong evidence for 
this. 

A summary of our results at the different dimensions 
are outlined in Tab. [H For the continuous phase tran- 
sitions, a generally accepted approach of obtaining more 
accurate values for the critical exponents is by examin- 
ing the local slopes (see for example [9]). For such an 
approach we plot, for example, S{t) against 1/t where 



Sit) 



ln[P{t)/P{t/m)] 
ln(m) 



(5) 



and 771 is the local range over which the slope is mea- 
sured. Due to the previously mentioned large corrections 
to scaling we observed, such an approach did not yield 
accurate results. For the method to work properly, the 
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sc(;2,iVoc,K) 


Order of 
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V 






phase transition 




(if applicable) (if applicable) 


1 




continuous 


0.071754 


0.160 


0.313 


log(32)/log(9) 


~ 1.5573 5'C(9^32,4) 


continuous 


0.08553 


0.30 


0.29 


log(16)/log(5) 


~ 1.7227 SC(5^16,6) 


continuous 


0.08819 


0.39 


0.24 


log(12)/log(4) 


~ 1.7925 5'C(4^12,6) 


? 


? 


? 


? 


log(8)/log(3) 


- 1.8928 S'C(3^ 8, 3/4/5) 


first-order 


0.093 


n/a 


n/a 


2 




first-order 


0.0973 


n/a 


n/a 



TABLE L Summary of the results. The data for d = 1,2 was obtained in Q]. Due to the large corrections to scaling, we were 
unable to provide strong evidence for the order of the phase transition for df ~ 1.7925. While we predict it to be continuous, 
we leave this row with question marks due to our uncertainty. 



number of time steps over which the data was collected 
would have to be raised significantly, increasing the re- 
quired computer time further. The values for the ex- 
ponents that we obtained then were estimated only by 
measurement of the gradient of the above plots and are 
therefore accurate to only ±0.005. 

A check for the consistency of the exponents can be ob- 
tained by the scaling relations. From the scaling relation 

i 

7^ + 26 = d/z, (6) 

we can check our values of rj and 5 using the fractal 
dimension for d and obtaining the dynamical exponent 
z = v\\lv^. This exponent can be obtained by examining 
the mean square distance of spread from the initial seed 
of the population R^{t) averaged over surviving runs 
only. At the critical point, we expect the asymptotic 
power law behavior R'^{t) 6c i^^^ Q. For df = 1.5573 we 
show the results in the inset of Fig. [3] a) where we plot 
R^{t) along with i^/^ using the value of z obtained from 
the scaling relation (O given by the obtained values 
for rj and S. We see very good agreement between the 
exponents. For larger dimensions, however, the crossover 



effects render, in particular R^{t), inaccurate. Oscil- 
lations in R^{t), perhaps due to the fractal structure, 
delay the power-law behaviour to large values of t that 
are impractical to simulate. 

Remarks: Having examined our model in different 
fractal dimensions, we found that the dimension at 
which the order of the phase transition changes is in 
the range 1.7227 < d{^ < 1.8928. We predict the lower 
bound of this range to be more accurately given by 
1.7927 but offer no firm proof. We predict that below 
df^ the larger correlations between the particles ensure 
that a continuous phase transition is observed. 

This present study was limited by the available com- 
puter power (over two years of processing time was re- 
quired in total). Due to the large corrections to scaling, 
it was difficult to obtain truly accurate results for the 
critical points and exponents in the continuous regime. 
Given more computer time, we would be able to run the 
simulations over a greater number of time steps to in- 
crease the accuracy. Further, it would also be interesting 
to examine if the values of the critical points and expo- 
nents depend on the type of fractal lattices used. 
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